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DE algorithm improves the DE by encoding individual parameters to produce 
and propagate better solutions. This paper proposes a self-adaptive differen- 
tial evolution algorithm with dynamic fitness-ranking mutation and pheromone 
strategy (SDE-FMP). The algorithm introduces the dynamical mutation opera- 
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Mutation strategy responding trial vector. The values are changed according to the pheromone 
Pheromone strategy when the trial vector is inferior in the selection, whereas the pheromone is in- 
Self-adaptation creased when the trial vector is superior. In addition, the algorithm also employs 


the resetting operation to unlearn and relearn the dominant pheromone values 
in the progressing search. The proposed SDE-FMP algorithm using the suitable 
resetting periods is compared with the well-known adaptive DE algorithms on 
several test problems. The results show that SDE-FMP can give high-precision 
solutions and outperforms the compared methods. 
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1. INTRODUCTION 

Solving multimodal, high-dimensional, and non-linear real-world optimization problems requires 
well-designed efficiency optimization methods [1]-[5]. To address this challenge, many researchers have pro- 
posed evolutionary algorithms such as genetic algorithm (GA) [6], ant colony optimization (ACO) [7], particle 
swarm optimization (PSO) [8], [Q], artificial bee colony algorithm (ABC) [T0], and differential evolution (DE) 
(1), for these problems. 

DE is a population-based global search algorithm introduced by Storn and Price in 1997 for continuous 
optimization. Its operations consist of mutation, crossover, and selection [11]. The performance of DE depends 
on the control parameters: scaling factor F and crossover rate CR. The scaling factor controls the step size 
of the mutation operation, and the crossover rate indicates the probability of exchanging elements between the 
mutant and target vectors. These control parameter values significantly affect the algorithm’s performance, and 
DE with the fixed parameters F and C'R are only suitable for specific problems. Thus, many control parameter 
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adaptation techniques have been proposed to improve the algorithm’s performance. 

The adaptive parameter technique uses the overall feedback from the search to adjust the parameter 
values, and the self-adaptive approach encodes the parameter values to the individuals and propagates them to 
others. These adaptation techniques can also modify the mutation operation. The DE algorithms with adap- 
tive parameters or improved mutation operations can accelerate the search process. However, they may face 
challenges in achieving high-precision solutions, such as getting stuck in local optimum or losing population 
diversity. Therefore, further improvements and strategies are necessary to enhance the adaptive approaches. 

This paper proposes a self-adaptive differential evolution algorithm with dynamic fitness-ranking mu- 
tation and pheromone strategy called SDE-FMP. The algorithm introduces a fitness-ranking mutation strategy 
that dynamically divides the population into three groups according to fitness rank. Then, it selects groups 
and their vectors with adaptive probabilities to create a mutant vector. The SDE-FMP also employs the self- 
adaptive control parameter adaptation by encoding the pre-assigned F and C'R values to the target vector and 
adjusting these values with the pheromone information. This self-adaptive process learns to find the appropriate 
parameters naturally. In addition, the algorithm incorporates the resetting operation to manage the dominant 
pheromone and enhance the efficiency of probabilities. The main contributions of our proposed algorithm are 
the ability to achieve high-precision solutions and superior performance to the compared methods. 

The remainder of this paper is organized as follows: i) section 2 briefly reviews some related work; 
ii) section 3 gives details of the proposed SDE-FMP algorithm; iii) the experiment sets are given in section 4; 
iv) section 5 discusses experimental results; v) the results and discussion are presented in section 6; and vi) the 
conclusion is given in section 7. 


2. LITERATURE REVIEW 

This section reviews the basic DE algorithm, self-adaptive DE algorithms, adaptive DE algorithms, 
and continuous ACO algorithms that inspire our proposed algorithm. 

The basic DE algorithm has three iterative operations: mutation, crossover, and selection. The al- 
gorithm generates NP population vectors from feasible solution space and indicates the best vector. For 
j = 1,..., NP, the mutation operation creates the mutant vector v; by randomly choosing three distinct 
vectors different from the target vector x; as (1): 


VU; = Er, + F- (Era — Erg) (1) 


where F is the scaling factor. 


Next, the crossover operation constructs the trial vector u; by exchanging the component of x; and vj 
with the crossover rate C'R. Finally, the selection operation updates the target vector with the trial vector when 
the fitness value of u; is better than that of x;. 

Since the basic mutation equation with the fixed parameters F and CR cannot solve a wide range 
of problems, the parameter adaptation and enhanced mutation strategy have been proposed to improve the 
performance of basic DE. 


2.1. Self-adaptive DE algorithms 

The self-adaptive DE algorithms encode the control parameter values to each target vector, adapt them 
based on the feedback of the search, and propagate the better ones to the next generation. 

The jDE by Brest et al. [13] is the self-adaptive DE algorithm encoding F; and C'R; to ith target 
vector with initial values 0.5 and 0.9, respectively. Mutation and crossover operations use these values from the 
target vector with the probabilities of 0.9; otherwise, it generates anew from [0.5, 0.9] and [0, 1], respectively. 
The algorithm does not use feedback from the selection operation. The jDE outperforms the basic DE with 
F = 0.5 and CR = 0.9 on several benchmark functions. Cheng et al. proposed DE with FDDE strategy 
that uses the combined fitness and diversity rankings to position the random vectors in a mutation operator 
where the diversity ranking computes from the difference of the median fitness and individual fitness values. 
The strategy improves the performance of jDE in both low-dimensional and high-dimensional problems. This 
work indicates the role of position selection for vectors in the mutation operator. Qin et al. proposed 
a differential evolution algorithm with strategy adaptation called SaDE. The algorithm uses four mutation 
operations and encodes their indices to each target vector. It assigns corresponding parameters F and C R from 
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Normal distributions and changes encoded information based on the probability of success and failure in the 
selection operation. The experimental results show that SaDE outperforms basic DE and jDE. The EPSDE 
algorithm by Mallipeddia et al. uses an ensemble of parameters and mutation strategies. This algorithm 
initially encodes a mutation strategy for a target vector and assigns corresponding parameters F and C'R. If the 
generated trial vector is not better than its target vector, the algorithm changes the target’s encoded information 
to the new one. EPSDE outperforms the classic DE and three adaptive DE algorithms: jDE, SaDE, and JADE. 


2.2. Adaptive DE algorithms 

The adaptive DE algorithm adjusts the control parameter with the overall feedback from the search. 
Then, the newly generated parameter will be biased according to the search feedback. 

Zhang and Sanderson presented the JADE algorithm that generates parameters F and CR from 
the Normal and Cauchy distributions, respectively. The algorithm implements an external archive for keeping 
the inferior vectors from selection and uses some top best vectors and the archived vectors in mutation. The 
results show that JADE is better than the classic DE and some adaptive DE algorithms. Tanabe and Fuku- 
naga introduced success-history-based parameter adaptation for differential evolution (SHADE), which 
improves the JADE algorithm with historical memories for each individual to update the mean of distributions. 
Experimental results show that SHADE is competitive with EPSDE, JADE, and CoDE. Next, Wang et al. 
introduced DE with composite trial vector generation strategies and control parameters called CoDE. It gen- 
erates three trial vectors with three different mutation operations and chooses the best one to compete with 
the target vector. The results show that CoDE is better than jDE, JADE, SaDE, and EPSDE algorithms. Zou 
et al. presented the CUSDE algorithm that uses a new mutation strategy by selecting the vectors with 
the probability calculated from the number of consecutive unsuccessful updates and removing individuals with 
those large numbers. The basic DE with this approach outperforms basic DE and some adaptive DEs. 


2.3. Ant colony optimization 

ACO is the population-based algorithm using the pheromone strategy to guide the ant population to 
locate optimal solutions for discrete optimization [21]. The ACO for continuous optimization requires dividing 
the initial spaces into discrete subspaces for formulating the pheromone structure. During the search process, 
the algorithm constructs a new solution vector with the components corresponding to subspaces according to 
the pheromone gathered from better solutions [22]. 

The ACOp algorithm is the first ACO algorithm for continuous optimization introduced by Dorigo 
and Socha [23]. The algorithm uses a pheromone structure to store the best solution components in an archive 
table and generates a new candidate solution with the corresponding distributions. At the end of the generation, 
it updates the pheromone by adding better candidate solutions and removing the worst archive solutions in the 
table. The performance of ACOR is competitive with other probability learning methods. 

Xiao and Li presented the hybrid of ACOpr and DE algorithms called HACO. It uses DE to 
generate new candidate solutions for the ACOR algorithm. The experimental results show that HACO per- 
forms better than ACO pg algorithms. Singsathid and Wetweerapong introduced the ACO with the domain 
partitioning technique called PACO. It generates solution components from partition points of adaptive search 
subspace that cover the best solutions and updates the pheromone according to the newly obtained best solution. 
The experiments show that PACO outperforms some well-known continuous ACO algorithms. 


3. THE PROPOSED SDE-FMP ALGORITHM 

We propose a self-adaptive differential evolution algorithm with dynamic fitness-ranking mutation 
and pheromone strategy, called SDE-FMP, which improve classic DE mutation by using the dynamic mutation 
and self-adaptive control parameters controlled with pheromone strategy and resetting operation. The details 
of SDE-FMP are as follows. 


3.1. New dynamic mutation strategy for SDE-FMP 

SDE-FMP sorts the population vectors at the beginning of each generation according to the fitness 
values and divides them into three groups with the same size, denoted by G1, G2, G3 where they represent the 
top best individuals, intermediate individuals, and worst individuals, respectively. 

To create a mutant vector for each target vector 7;,2 = 1,2,..., NP, the algorithm chooses a group 
gk (from G1, G2, or G3) for each position k = 1, 2,3 of the mutation equation (g2 and g3 must be different 
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for diversity) with probability vectors PropR,, calculated from pheromone vectors Pheromone Rpg. Note that 
SDE-FMP initializes all components of pheromone vectors to be 1, ensuring equal probability. 

Next, the algorithm chooses three distinct random vectors £R, R3, £R, Which differ from x; by 
uniformly selecting random vectors in the corresponding selected g% groups to generate a mutant vector as (2): 


vi = ZR, + F(t): (tr, — ERs) 4 


where the scaling factor F (i) corresponds to target vector xi. 
The mutant vector enters the crossover operation to generate the trial vector u; as (3): 


ij isj ; j= 
uij = Vig ' SJ CR(i) or J Irana GB) 
Zij 3 otherwise 
where j = 1,..., D; sj is a uniform random number in (0, 1) and Iyana is a randomly fixed integer from 1 to 


D. 

At the selection operation, if u; is better than x;, the algorithm updates the pheromone at the lth 
position of gg = G by adding one to the associated pheromone vector position to reinforce the pheromone 
information related to a successful solution as (4): 


PheromoneR;(l) = PheromoneR, (1) +1 (4) 


At the end of each generation, the algorithm normalizes the Pheromone Rygę to the probability vector 
PropRyx as (5): 


Ph l 
Prok E eromoneR, (l) 


5 
Xn- PheromoneR,(m) a 


for l = 1,2,3. 


3.2. Self-adaptive control parameters of F and CR 

We use scaling factor values F = 0.5,0.7,0.9 to control the step size for the mutant vectors and 
crossover rate values CR = 0.1, 0.9 to balance between intensifying and diversifying the search. CR = 0.1 
is suitable for a local search, while C'R = 0.9 is suitable for a global search. So, we have six combinations of 
these values. 

At initialization, the algorithm sets all components of the pheromone vector Pheromone FC R to be 
one and encodes a random pair of F (i) and CR(i) for each target vector x;. The algorithm uses F (i) and 
C R(i) from z; to generate a trial vector u; from mutation and crossover operations. 

At the selection operation, if u; is better than z;, the target vector retains its current F' (i) and CR(i) 
values and the associated pheromone vector value Pheromone FC R(t) is incremented by one. 


PheromoneF'CR(t) = PheromoneF'CR(t) +1 (6) 


where ¢ is the corresponding index of that combination. Otherwise, the target vector is re-encoded with a 
new random pair of F' (i) and CR(i) based on the probability vector PropF C R, allowing the target vector to 
explore new parameter values. Note that PropF'CR is calculated from Pheromone FC R for each generation 
as (7): 


PheromoneF'CR(t) 
Si _, PheromoneFCR(m) 


m=1 


PropFCR(t) = (7) 


fort = 1,2,...,6. 


3.3. The pheromone resetting 

At the end of each generation, SDE-FMP determines whether the pheromone vectors have reached the 
specified thresholds to prevent the dominance of certain pheromone values and promote fair competition among 
choices. The two parameters rg and rp represent the thresholds for PheromoneR, and PheromoneF'CR, 
respectively. If the sum of any pheromone vector is greater than the corresponding threshold, the algorithm 


Bulletin of Electr Eng & Inf, Vol. 13, No. 1, February 2024: 559-571 


Bulletin of Electr Eng & Inf ISSN: 2302-9285 o 563 


resets all elements of those pheromone vectors back to 1. The pseudo code of SDE-FMP is presented in 
Algorithm 1. 


Algorithm 1 SDE-FMP algorithm 


1: Initialize the population of N P individuals 
2: Find the best vector xyes; and its best function value frest 
3: Encode the F (i), CR(i) values to each target vector z;, i = 1,..., NP 
4: Set all elements of PheromoneRy,, k = 1,2,3 and PheromoneF'CR to 1 
5: Set all elements of ProbR;, to be equal for k = 1, 2,3 
6: Set all elements of ProbF' CR to be equal 
7: Set number of function evaluations n f = 0 
8: Set the VTR or the maximum number of function evaluations marnf 
9: while stopping condition is not satisfied do 
10: Sort the population individuals according to fitness ranking and divide them into G1, G2, G3 groups 
11: for i = 1 : NP do 
12: Choose distinct xp, vectors with the probabilities ProbRg, k = 1,2,3 
13: Generate a mutant vector using eq. 
14: Apply the crossover operation eq. to get a trial vector u; 
15: Evaluate f(u,;) andnf nf +1 
16: if f(u) < f(a,) then 
17: Replace x; with u; 
18: Update each PheromoneRy, using eq. 
19: Update PheromoneF'CR using eq. (6) 
20: if f(ui) < foes: then 
21: Replace £pest with u; 
22: end if 
23: else 
24: Re-encode the new random F (i), CR(i) for x; according to ProbF CR 
25: end if 
26: end for 
27: if X PheromoneR, > rg for some k then 
28: PheromoneR, + 1, k = 1,2,3 
29: end if 
30: if $, PheromoneFCR > r, then 
31: PheromoneFCR + 1 
32: end if 


33: Normalize Pheromone Rx to be ProbR, using eq.(5) 

34: Normalize PheromoneF'CR to be ProbF'CR using eq.(7) 
35: end while 

36: Report the obtained pest, fbest and nf 


4. EXPERIMENTAL DESIGN 

The performance of SDE-FMP is tested on eight selected benchmark functions that cover four main 
types: uni-modal, multi-modal, separable, and non-separable. Their formulae and search ranges are presented 
in Table 1. First, we design a preliminary experiment for finding the suitable values rg and rp for all types of 
problems. Then, we conduct two comparison experiments to evaluate the performance of SDE-FMP against 
other adaptive DE algorithms. The details of each experiment are given in the following subsections. 
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Table 1. Test functions 


Function Formula Search range 
D 
Sphere Fi (z) = X (ey? —100, 100 
i=1 
D i 
Schwefel 1.2 F2(x) = XCO z)? —100, 100 
i=1 j=1 
D-1 
Rosenbrock F3(x) = (100(a441 — £2)? + (xi — 1)”) —100, 100 
i=1 
i 1 2 D a 
Griewank Falz) = x? cos(—=) +1 —600, 600 
o- ag 8 - Tom 
D 
Rastrigin Fs5(x) = X `[z? — 10cos(2ra;) + 10] —5.12, 5.12] 
i=1 
12 
Ackl F6(x) = —20e 0.2 2 —32.32 
ckley (a exp( 5) 


D 


1 
— exp(5 a cos(2rx;)) +20 +e 
i=1 


D 
Schwefel F;(x) = 418.98288727243369 - D — X` (x; sin(vlæ:])) — [—500, 500] 
i=l 
D D i 
Schwefel 2.22 Fg(x) = X` |as| + [| lz: [—100, 100] 
i=1 i=l 


4.1. Finding the suitable values r, and r, for SDE-FMP 

To find the suitable values for resetting periods r4 and rp, the dimensions of test functions are D = 
10, 30. The number of population NP = 30, magn f = 20000D and VTR = 10~!° are used. The parameters 
rg and rp are varied as rọ = 200, 500 and r, = 200, 300, 400. The algorithm performs 50 independent runs for 
each configuration. We report the number of successful runs (NS), the mean number of function evaluations 
(meanNF), and the percentage of the standard deviation of function evaluations (%SD). 


4.2. Comparing the performance of SDE-FMP with other adaptive DE algorithms using VT R 

We use the obtained values rg and rp to compare the performance of SDE-FMP with some well-known 
adaptive DE algorithms: JADE [17], CoDE [19], jDE [13], and SaDE [15]. The SDE-FMP uses the same 
setting as the first experiment, while the compared algorithms use the parameter settings as in their original 
papers. The dimensions are varied as D = 10,30,50. All algorithms perform 100 independent runs for 
each configuration. The MATLAB source codes of JADE, CoDE, jDE, and SaDE are available from Zhang’s 
homepage: http://dces.essex.ac.uk/staff/qzhang. 


4.3. Comparing the performance of SDE-FMP with other adaptive DE algorithms using maxnf on CEC 
2005 benchmark functions 

We compare the performance of SDE-FMP using suitable r, and r, with SaDE [15], FIDE F [14], 
and CUSDE on 30-dimensional benchmark functions of CEC 2005 [26]. The experiment reports the 
mean of optimal values and standard deviation using mazn f = 10000D over 50 independent runs. The best 
values of the compared algorithms are from their original papers. We use the t-test at the significance level of 
0.05 to compare their performances. The symbols +, 0, — represent that the mean of the optimal value of the 
SDE-FMP is superior to, equal to, and inferior to the compared algorithm, respectively. 


5. EXPERIMENTAL RESULTS 
5.1. The suitable values r, and r, for SDE-FMP 

The experiment finds the suitable values 7, and rp that give the highest number of successful runs 
and the lowest meanNF for SDE-FMP. Table 2 shows that three combinations of (rg, rp) i.e., (200, 400), (500, 
200), and (500, 300) give 50 successful runs for all cases. We highlight the lowest meanNF values among these 
three combinations for each case and obtain r = 500 and r, = 300 as the best values. 
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Table 2. The performance comparison of SDE-FMP with different values r4 and r, over 50 independent runs 


i Tg =200 Tg =500 Tg = 1000 
as Statistics rp =200 rp=300 rp=400 rp=200 rp=300 rp=400 rp=200 rp=300 rp = 400 
NS 50 50 50 30 50 50 50 50 50 

10 meanNF 9122.62 8966.98 8939.94 8589.55 8569.9 8593.13 8309.1 8301.7 8330.37 

n %SD 5.54 5.65 6.66 8.29 7.00 6.32 5.64 9.38 8.78 
NS 50 50 50 50 50 50 50 50 50 

30 meanNF 31727.14  32191.1  32070.8 31035.48 30801.56 31003.33 26165.16 29451.24 29621.88 

%SD 4.24 3.08 3.17 3.87 4.21 4.05 5.54 5.65 6.13 

NS 50 50 50 50 50 50 50 50 50 

10 meanNF 148864 148136 15007.06 14532.58 14427.32 14364.64 13775.38 14596.9 14374.22 

ji %SD 6.98 8.68 9.23 11.01 10.5 11.33 7.50 10.68 12.54 
NS 50 50 50 50 50 50 50 50 50 

30 meanNF 141685.6 139038.04  139334.24 133898.42 131200.62 129951.46 107370.16 129738.66  127754.88 

%SD 6.69 8.28 6.92 8.00 9.47 5.56 12 8.62 10.38 

NS 48 47 50 50 50 46 44 50 48 

10 meanNF 25334.52 25313.85 25288.82 24499,9  24470.8 24742.02 216139 24017.18 23592.12 

P %SD 9.77 13.12 10.07 9.77 14.57 11.14 10.03 14.55 17.65 
3 NS 49 49 50 50 50 50 50 47 47 
30 meanNF 219922.79 227484.93 246704.98 202595.87 207844.62 218021.04 144629.44 192881.38 199274.53 

%SD 103 113 11.51 9.87 12.76 10.95 8.48 12.23 12.35 

NS 50 50 50 50 50 50 50 49 49 

10 meanNF 23443.48 24055.58 23895.28 22821 23882.22 24615.22 23696.84 23748.12  23536.1 

s %SD 17.16 14.51 18.72 19.98 20.50 24.70 19.75 16.59 17.67 
NS 50 50 50 50 50 49 50 50 50 

30 meanNF 35335.16  34497.8 34799.02 33826.88 33284.62 34164.87 28456.36 32186.22  32371.8 

%SD 5.73 5.70 7.00 6.83 5.98 7.32 8.33 6.47 8.56 

NS 50 50 50 50 50 50 50 50 50 

10 meanNF 12863.76  12887.7 12956.24  12571.5 12342.33 12494.54 12598.32 12326.33 12362.02 

m WSD 4.95 4.26 4.36 4.96 5.58 5.29 5.35 6.94 6.32 
; NS 50 50 50 50 50 50 50 50 50 
30 meanNF 58730.32 57263.44 57201.62 57217.98 56480.04 55874.78 52726.84 55993.94 55891.36 

%SD 4.01 3.22 3.54 4.15 3.92 4.37 5.36 4.18 3.72 

NS 50 50 50 50 50 50 50 50 50 

10 meanNF 15230.08 15207.54 5109.3 14476.26 14143.74 14008.16 13817.44 13791.78 13784.34 

m WSD 4.99 5.05 5.38 5.35 4.73 5.63 7.42 8.34 8.84 
NS 50 50 50 50 50 50 50 50 50 

30 meanNF 520431  51920.1 51532.06 50354.46  49440.6 50054.94 48582.88 48513.56 48565.88 

%SD 2.30 2.44 2.36 2.92 4.47 2.54 4.28 5.06 4.92 

NS 50 50 50 50 50 50 50 50 50 

10 meanNF 11398.32 11515.02 11405.52 11092.34 10859.18 10954.26 10957.84 10834.18 10900.8 

= %SD 4.33 4.89 4.87 5.48 5.72 6.36 5.27 6.55 6.51 
NS 50 50 50 50 50 50 50 48 50 

30 meanNF 41162.04 41313.46  40652.2 40615.06 39669.68 39620.16 35394.56 39113.77 38681.88 

%SD 2.88 2.75 2.79 3.74 4.63 4.70 6.57 5.46 5.39 

NS 50 50 50 50 50 50 50 50 50 

10 meanNF 15752.34 15797.36 15921.62 15143.78 15199.58 15148.26 14683.66 14625.28 14772.74 

r %SD 4.09 4.03 4.34 5.04 4.33 5.93 6.75 5.50 6.02 
NS 50 50 50 50 50 50 50 50 50 

30 meanNF 53337.04 53095.36 52594.54 51815.44 51384.32 51184.22  50815.5 50257.46 50110.68 

%SD 221 2.28 1.63 2.45 2.69 2.53 3.40 3.16 3.62 

NO. 50 successful runs cases 14/16 14/16 16/16 16/16 16/16 14/16 15/16 13/16 13/16 


5.2. Performance comparison of SDE-FMP with other adaptive DE algorithms using VTR 

We compare the performance of SDE-FMP with SaDE, CoDE, jDE, and JADE using the VTR = 
10710, Table 3 presents the number of successful runs and meanNF and highlights the best results that give 
100 successful runs with the lowest meanNF. 

The results show that SDE-FMP, SaDE, CoDE, jDE, and JADE give 100 successful runs for 24, 12, 16, 
21, and 17 cases, respectively. Their lowest mean counts are 18, 0, 1, 0, and 5, respectively. Therefore, SDE- 
FMP achieves the best performance. Note that all compared algorithms cannot achieve high-quality solutions 
for the Rosenbrock function F3. 
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Table 3. The performance comparison of SDE-FMP, SaDE, CoDE, jDE, and JADE over 100 independent runs 

F D Statistics SDE-FMP SaDE CoDE jDE JADE 

10 NS 100 100 100 100 100 

meanNF(%SD) 8549(7.70) 12304(2.82)  24199(2.80) 28358(2.57)  23035(2.95) 

F, 30 NS 100 100 100 100 100 

meanNF(%SD) 30732(4.78) 31005(3.94)  61089(2.32) 67349(1.37)  34577(3.07) 

50 NS 100 100 100 100 100 

meanNF(%SD) 54935(2.62) 56234(3.83)  87384(2.33) 95734(1.87)  42939(2.52) 

10 NS 100 100 100 100 100 

meanNF(%SD) 14372(9.92)  26321(19.18)  38828(4.08) 65857(4.57)  29031(5.31) 

F> 30 NS 100 100 100 100 100 

meanNF(%SD) 131979(7.96)  394091(6.68) | 208775(1.00)  377665(5.12)  87170(4.94) 

50 NS 100 0 9 22 100 

meanNF(%SD)  451807(16.85) 0(0.00) 486460(2.19)  970940(2.44) 218639(4.31) 

10 NS 100 38 100 100 T 

meanNF(%SD) 24218(13.20) 14223(25.84) 58911(5.01) 137735(12.69) — 53357(12.05) 

Fs 30 NS 100 0 81 31 10 

meanNF(%SD) 208507(11.35) 0(0) = 281516(6.59) 529180(4.71) — 126790(3.09) 

50 NS 100 0 1 4 4 

meanNF(%SD) 769397(11.69) 0(0) 488460(0) 946775(6.24) 225125(3.2) 

10 NS 100 100 100 100 100 

meanNF(%SD) 23283(23.00) 23816(1.84)  13619(1.87) 47520(1.72)  54266(3.47) 

F, 30 NS 100 100 100 100 100 

meanNF(%SD) 33712(6.66) 53759(2.44) 122277(1.81) 114079(1.52) 198059(1.48) 

50 NS 100 100 100 100 100 

meanNF(%SD) 56498(3.67) 87137(3.90) 174414(2.19)  165118(1.52)  314970(7.85) 

10 NS 100 100 100 100 100 

meanNF(%SD) 12460(5.59) 24755(4.07)  36493(3.29) 44271(3.47)  47759(1.95) 

Fs 30 NS 100 97 97 100 100 

5 meanNF(%SD) 56188(4.13) 73422(5.32) 150794(4.22)  123334(2.68) 143686(1.30) 

50 NS 100 57 49 100 100 

meanNF(%SD) 148488(4.75)  133396(5.165) 290445(5.28)  193196(2.80)  236536(1.4) 

10 NS 100 100 100 100 49 

meanNF(%SD) 14197(6.54) 20601(3.93)  28871(3.29) 32980(2.81)  43536(4.79) 

Fg 30 NS 100 100 99 100 55 

meanNF(%SD) 49699(3.46) 53811(3.82)  92755(3.56) 84608(2.64)  138990(1.71) 

50 NS 100 91 100 100 16 

meanNF(%SD) 87009(2.22) 97754(3.76) 163292(4.12)  126238(2.29) 227293(2.41) 

10 NS 100 100 100 100 100 

meanNF(%SD) 11002(6.84) 20237(2.00) 40587(2.00) 47096(1.50)  40421(3.12) 

F; 30 NS 100 89 100 100 100 

meanNF(%SD) 40029(3.77) 68850(1.32) 98994(1.90) 106889(1.67) 55491(2.99) 

50 NS 100 12 100 100 100 

meanNF(%SD) 76032(2.61)  324766(1.08) 138260(1.75)  149720(1.41) —-67224(2.06) 

10 NS 100 100 100 100 100 

meanNF(%SD) 15243(5.22)  36324(15.63)  59423(6.46)  65944(10.15)  86969(5.53) 

Fe 30 NS 100 67 100 100 100 

meanNF(%SD) 51775(2.66) 32029(5.01)  65511(7.48) 70432(3.23)  41502(28.77) 

50 NS 100 52 95 100 8 

meanNF(%SD) 89138(1.60) 56207(3.84)  88634(3.51) 97001(2.43)  43912(2.36) 

NO. 100 successful runs cases 24/24 12/24 16/24 21/24 17/24 

NO. lowest meanNF cases 18/24 0/24 1/24 0/24 5/24 


5.3. Performance comparison of SDE-FMP with other adaptive DE algorithms using maxnf on the CEC 
2005 benchmark functions 

We compare the performance of SDE-FMP with SaDE, FDDE_F, and CUSDE using the mean of 
obtained best values. Table 4 shows that the superior cases of SDE-FMP to SaDE, FDDE_F, and CUSDE are 
12, 9, and 10, whereas the inferior ones are 5, 6, and 5. Therefore, SDE-FMP overall outperforms the compared 
methods on CEC 2005 benchmark functions. 
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Table 4. The performance comparison of SDE-FMP, SADE, FDDE_F, and CUSDE on 30-dimensional 
CEC2005 benchmark functions 


F SDE-FMP SaDE FDDE_F CUSDE 
MEAN std. MEAN std. MEAN std. MEAN std. 
fl 0.00E+00 0.00E+00 0.00E+00 = 0.00E+00 1.40E-03 +  1.47E-03 6.73E30 + 5.62E-29 
f2 1.48E-14 2.52F-14 2.77E-06 + 8.52E-06 5.20E-13 +  3.36E-13 2.82E-08 + 2.28E-08 
f3 3.58E+04 3.07E+04 5.33E+05 + 4.34E+05 1.01E+05 + 8.22E+04 2.32E+05 + 1.09E+05 
f4  9.77E-01 6.90E+00 1.93E+02 + 3.22E+02 2.00E-03 - 1.63E-03 1.32BE-05 -  1.10E-05 
f5 6.83E+02 4.43E+02 3.76E+03 + 6.12E+02 1.80E+02 - 1.47E+02 9.05E-05 -  7.68E-05 
f6 1.75E-13 6.75E-13  5.28E+01 + 4.15E+01 6.22E-01 + 5.07E+00 5.47E-09 + 2.60E-08 
f7  6.84B-03 9.19F-03  1.65E-02 + 1.58E-02 7.20E-02 + 5.87E-02 0.00E+00 - 0.00E+00 
f8 2.09E+01 4.09F-02 2.09E+01 =  1.58E-02 2.09E+01 = 1.70E-02 2.09E+01 = 5.53E-02 
f9  0.00E+00 0.00E+00 0.00E+00 = 0.00E+00 5.45E+01 + 4.44E+00 5.23E+01 + 2.13E+01 
flO 1.08E+02 1.55E+01 4.76E+01 - 1.26E+01 4.66E+01 - 3.80E+00 1.72E+02 - 8.29E+00 
fll 2.92E+01 1.32E+00 1.68E+01 - 1.64E+00 2.79E+01 = 2.27E+00 3.73E+01 +  7.16E-01 
f12 3.43E+04 5.37E+03 3.44E+03 - 4.42E+03 3.74E+03 - 5.05E+03 3.59E+03  - 4.99E+03 
f13  2.83E+00 2.14E-01 3.84E+00 + 2.66E-01 1.67E+00 - 1.36E+00 141E+01 + 847E-01 
fl4 1.29E+01 1.74E-01 1.26E+01 = 2.83E-01 1.29F+01 = 1.05E+00 1.30E+01 = 2.12E-01 
f15 2.10E+02 2.67E+01 3.85E+02 + 4.42E+01 4.00E+02 + 3.26B+01 3.33E+02 +  41.11E+02 
fl6 1.57E+02 5.81E+01 8.65E+01 - 5.65E+01 1.13E+02 + 9.21E+00 2.37E+02 + 8.45E+01 
f17 1.88E+02 5.11E+01 8.15E+01 - 3.46E+01 2.64E+02 + 2.15E+01 2.36E+02 + 5.98E+01 
f18 9.05E+02 1.53E+01 8.73E+02 = 5.44E+01 9.04E+02 = 2.37E+00 9.03E+02 =  1.18E-01 
f19 9.03E+02 2.13E+01 8.74E+02 = 5.44E+01 9.05E+02 = 7.37E+00 9.03E+02 =  1.89E-01 
f20 9.07E+02 1.35E+00 8.81E+02 + 5.22E+01 5.54E+02 - 4.54B-01 9.03E+02 =  7.99E-01 
f21 5.00E+02 1.83E-13 5.45E+02 + 2.15E+02 5.00E+02 = 0.00E+00 5.00E+02 =  1.62E-13 
f22 9.12E+02 1.02E+01 9.21E+02 + 2.66E+01 8.72E+02 = 7.10E-01 8.61E+02 = 3.77E+00 
f23 5.34E+02 3.98E-04 5.34E+02 = 8.27E-04 5.34E+02 = 4.35E-02 5.34E+02 = 3.94B-04 
f24 2.00E+02 8.66E-13 2.00E+02 = 8.54E-14 2.00E+02 = 0.00E+00 2.00E+02 = 2.94E-14 
f25  2.11E+02 7.32E-01 2.14E+02 + 2.35E+00 2.10E+02 = 1.72E-01 2.09E+02 = 2.46E-01 
(+/ = /-) 12/8/5 9/10/6 10/10/5 


6. DISCUSSION 

The SDE-FMP algorithm uses the pheromone strategy to adapt the probabilities for selecting sub- 
groups in mutation where a high pheromone indicates a suitable group for the corresponding position in the 
mutation equation. Then, the mutant vector has more potential to create a better trial vector during the search. 
The algorithm also uses the pheromone for self-adaptive control parameters F and C'R, where the most suc- 
cessful pair of F and CR has more potential to propagate to the next generations. The pheromone reset- 
ting is employed to eliminate the dominance and balance the cycle of gathering the pheromone and using the 
pheromone to improve the search performance. We obtain the suitable resetting periods rg = 500 and r, = 300 
for PheromoneR,, k = 1,2,3 and PheromoneF CR, respectively. 

We further investigate the impact of SDE-FMP’s features, including the dynamic mutation strategy, 
self-adaptive control parameters, and pheromone resetting. We compared the performance of the SDE-FMP 
with the SDE-FMP without the proposed mutation strategy (using basic DE mutation strategy), the SDE-FMP 
without self-adaptive control parameters (using fixed F = 0.5 and CR = 0.9 values), and the SDE-FMP 
without pheromone resetting on the 30-dimensional Rosenbrock function for 30 independent runs. The results 
presented in Table [5] demonstrate that the SDE-FMP significantly outperforms the SDE-FMP without each 
feature. The dynamic mutation strategy and self-adaptive control parameters play crucial roles in improving 
the convergence speed, while the pheromone resetting further enhances the achievement of high-precision 
solutions. 

Figures [IJand[2|illustrate the convergence graphs of SDE-FMP compared with SaDE, CoDE, jDE, and 
JADE on the Sphere, Griewank, Ackley, and Schwefel functions for 10 and 30 dimensions. They show that 
SDE-FMP can solve problems of various types faster than the compared algorithms. 


Self-adaptive differential evolution algorithm with dynamic fitness-ranking mutation ... (Pirapong Singsathid) 


568 o 


ISSN: 2302-9285 


Table 5. Comparison of SDE-FMP and SDE-FMP without dynamic mutation strategy, self-adaptive control 
parameter, and pheromone resetting on the 30-dimensional Rosenbrock function 


i SDE-FMP without SDE-FMP without SDE-FMP without 
Algorithm SDE-FMP ; : í : 
dynamic mutation strategy  self-adaptive control parameters pheromone resetting 
NS 30 4 8 23 
meanNF(%SD) 207681.97(9.52) 412309.50(35.91) 487927.25(20.31) 


257962.26(43.18) 
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Figure 1. Convergence graphs of SDE-FMP, SaDE, CoDE, jDE, and JADE on 10-dimensional sphere, 
Griewank, Ackley, and Schwefel functions 
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Figure 2. Convergence graphs of SDE-FMP, SaDE, CoDE, jDE, and JADE on 30-dimensional sphere, 
Griewank, Ackley and Schwefel functions 


7. CONCLUSION 

This paper presents a self-adaptive differential evolution algorithm with dynamic fitness-ranking mu- 
tation and pheromone strategy called SDE-FMP. The pheromone strategy manages the adaptive probabilities 
for a dynamic mutation and self-adaptive control parameters, and the resetting operation helps the search to 
prevent premature convergence and stagnation. As a result, the proposed algorithm can solve problems of var- 
ious types. Experiments indicate that SDE-FMP gives high-precision solutions and outperforms the compared 
methods on benchmark functions. 
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